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Abstract 

We present, theoretical predictions and Monte Carlo simulations, for a simple three 
matrix model that exhibits an exotic phase transition. The nature of the transition is 
very different if approached from the high or low temperature side. The high temperature 
phase is described by three self interacting random matrices with no background spacetime 
geometry. As the system cools there is a phase transition in which a classical two-sphere 
condenses to form the background geometry. The transition has an entropy jump or latent 
heat, yet the specific heat diverges as the transition is approached from low temperatures. 
We find no divergence or evidence of critical fluctuations when the transition is approached 
from the high temperature phase. At sufficiently low temperatures the system is described 
by small fluctuations, on a background classical two-sphere, of a U{1) gauge field coupled 
to a massive scalar field. The critical temperature is pushed upwards as the scalar field 
mass is increased. Once the geometrical phase is well established the specific heat takes 
the value 1 with the gauge and scalar fields each contributing 1/2. 

THE AUTHOR B. Y DEDICATES THIS WORK TO THE MEMORY OF HIS DAUGHTER 
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1 Introduction 



All the fundamental laws of physics are now understood in geometrical terms. The classical 
geometry that plays a fundamental role in our formulation of these laws has been vastly extended 
in noncommutative geometry [1]. However, we still have very little insight into the origins of 
spacetime geometry itself. This situation has been undergoing a significant evolution in recent 
years and it now seems possible to understand geometry as an emergent concept. The notion 
of geometry as an emergent concept is not new, see for example [2] for an inspiring discussion 
and [3, 4] for some recent ideas. We examine such a phenomenon in the context of a simple 
three matrix model [5-7]. 

Matrix models with a background noncommutative geometry have received attention as an 
alternative setting for the regularization of field theories [8-11] and as the configurations of DO 
branes in string theory [12,13]. In the model studied here, the situation is quite different. It has 
no background geometry in the high temperature phase and the geometry itself emerges as the 
system cools, much as a Bose condensate or superfluid emerges as a collective phenomenon at 
low temperatures. The simplicity of the model we study here allows for a detailed examination 
of such an exotic transition. We suspect the characteristic features of the transition may be 
generic to this novel phenomenon. 

In this article we study both theoretically and numerically a three matrix model with global 
5*0(3) symmetry whose energy functional or Euclidean action functional (see eq. (13.61) ) is a 
single trace of a quartic polynomial in the matrices Da- The model contains three parameters, 
the inverse temperature /? = = g'"^, and parameters m and /x which provide coefficients for 
the quartic and quadratic terms in the potential. 

We find that as the parameters are varied the model has a phase transition with two clearly 
distinct phases, one geometrical the other a matrix phase. Small fluctuations in the geometrical 
phase are those of a Yang-Mills and a scalar field around a ground state corresponding to a 
round two-sphere. In the matrix phase there is no background spacetime geometry and the 
fluctuations are those of the matrix entries around zero. In this article we focus on the subset 
of parameter space where in the large matrix limit the gauge group is Abelian. 

For finite but large N, at low temperature, the model exhibits fiuctuations around a fuzzy 
sphere [14]. In the infinite N limit the macroscopic geometry becomes classical. As the tem- 
perature is increased it undergoes a transition with latent heat so the entropy jumps, yet the 
model has critical fiuctuations and a divergent specific heat. As this critical coupling is ap- 
proached the fuzzy sphere radius expands to a critical radius and the sphere evaporates. The 
neighbourhood of the critical point exhibits all the standard symptoms of a continuous 2nd 
order transition, such as large scale fiuctuations, critical slowing down (of the Monte Carlo 
routine) and is characterized by a specific heat exponent which we argue is a = 1/2, a value 
consistent with our numerical simulations. In the high temperature (strong coupling) phase 
the model is a matrix model closely related to zero dimensional Yang-Mills theory. As the 
transition is approached from within this phase we find no evidence of critical fiuctuations and 
no divergence in the specific heat. 
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In the geometrical sphere phase the gauge couphng constant is — (5~^ and m parameterizes 
the mass of the scalar field. For small we find a transition with discontinuous internal 
energy U —< S > / (3 (so the entropy jumps across the transition [15]) while the specific heat 
is divergent as the transition is approached from the low temperature fuzzy sphere phase but 
finite when approached from the high temperature matrix phase. The fuzzy sphere emerges 
as the low temperature ground state which as one expects is the low entropy phase and the 
transition is characterized by both a latent heat and divergent fluctuations. 

To our knowledge it is the first clear example of a transition where the spacetime geometry 
is emergent. This transition itself is extremely unusual. We know of no other physical situation 
that has a transition with these features. Standard transitions are very dependent on the 
dimension of the background spacetime and when this is itself in transition an asymmetry of 
the approach to criticality is not so surprising. 

By studying the eigenvalues of operators in the theory we establish that, in the matrix 
phase, the matrices Da are characterised by continuous eigenvalue distributions which undergo 
a transition to a point spectrum characteristic of the fuzzy sphere phase as the temperature 
is lowered. The point spectrum is consistent with Da — L^/R where La are su{2) angular 
momentum generators in the irreducible representation given by the matrix size and R is the 
radius of the fuzzy sphere. The full model received an initial study in [7] while a simpler version, 
invariant under translations of arises naturally as the configuration of DO branes in the 
large k limit of a boundary Wess-Zumino-Novikov-Witten model [13] and has been studied 
numerically in [5]. If the mass parameters of the potential are related to the matrix size, the 
model becomes that introduced in [16]. The interpretation presented here is novel, as are the 
results on the entropy and critical behaviour and the extension to the full model. 

A short description of the results obtained in this article is given in [15]. This article is 
organised as follows. In section 2 we review the fuzzy sphere and its geometry. In section 3 we 
derive several theoretical predictions in the fuzzy sphere phase including the critical behaviour 
of the model. In section 4 we discuss the non-perturbative phase structure (phase diagram) 
and Monte Carlo numerical results for various observables. We conclude in section 5 with some 
discussion and speculations. 

2 The fuzzy sphere 

The ordinary round unit sphere, 5*^, can be defined as the two dimensional surface embedded 
in flat three dimensional space satisfying the equation X1h=i "a ~ ^ with nER^. One can 
use the nonholonomic coordinate system for the sphere. In this coordinate system 

a general function can be expanded as /(n) = Yl'^o fimyim{n) , where Yim are the standard 
spherical harmonics. The basic derivations are provided by the SO (3) generators deflned by 
Ca — —i^abciT'bdc and the Laplacian is — CaCa with eigenvalues l{l + 1), I — 0, ...,oo. 
Following Frohlich and Gaw§dzki [25] (or Connes [1] for spin geometry) the geometry of a 
Riemannian manifold can be encoded in a spectral triple. For the ordinary sphere the triple 
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is (C°°(S^), 7Y, where C°°(S'^) is the algebra of all functions f{n) on the sphere and H is 
the infinite dimensional Hilbert space of square intcgrable functions. Such a spectral triple is 
precisely the data that enters the scalar field action on the manifold. So a scalar field action 
which includes an appropriate Laplacian can specify the geometry. 

The fuzzy sphere can be viewed as a particular deformation of the above triple which is 
based on the fact that the sphere is the coadjoint orbit SU{2)/U{1) [14], 

gcTsg-' = riaaa , geSU{2) , neS\ (2.1) 

and is therefore a symplectic manifold which can be quantized in a canonical fashion by simply 
quantizing the volume form 

u — sin9d9/\d(l) — -eabc'^adnb/Kdrbc- (2-2) 

The result of this quantization is to replace the algebra C°°(5'^) by the algebra of NxN matrices 
Mat]\f. MatN becomes the A^^— dimensional Hilbert space structure when supplied with 
inner product {f,g) = j^Tr{f^g) where f,gEMatN. The spin IRR of su{2) has both a 
left and a right action on this this Hilbert space. For the left action the generators are La and 
satisfying [La, Lb] — iSabcLc, and X^a-^a = C2 — ^^^^ ■ The spherical harmonics yjr„(n) become 
the canonical su{2) polarization tensors Yi^ and form a basis for H^. These are defined by 

[La, [La, Yi^]] = l{l + l)Yi^ , [L±, Yi^] = ,/{lTm){l±m + l)Yi^±i , [L3, Yim] = mYim- (2.3) 

and satisfy 



and the completeness relation 



N-l I 



E E ^L^lrn = 1 (2.5) 
;=0 m=-l 

The "coordinates functions" on the fuzzy sphere 3% are defined to be proportional to Yi^ 
tensors (as in the continuum) and satisfy 

xj + xl + xl^l , [xa, Xb] = -^eabcXc, whcrc Xa = (2.6) 

V C2 V C2 



"Fuzzy" functions on S% are elements of the matrix algebra while derivations are inner and 
given by the generators of the adjoint action of su{2) defined by Ca(f) '■= [La,4']- A natural 
choice of the Laplacian on the fuzzy sphere is therefore given by the Casimir operator 

jC,'^[La,[La,..]]. (2.7) 

Thus the algebra of matrices Matj^ with = L + 1 decomposes under the action of su{2) as 
= 0©1©2©..©L, with the first ^ standing for the left action while the other ^ stands 



for the right action of sw(2). It is not difficult to convince ourselves that this Laplacain has a 
cut-off spectrum with eigenvalues /(/ + 1) where I = 0, 1, L. Given the above discussion we 
see that a general fuzzy function (or element of the algebra) on Sfj can be expanded in terms 
of polarization tensors as follows / = Xlz^o XlL=-« fimXim- The continuum limit is given by 
L — >oo. Therefore the fuzzy sphere can be described as a sequence of triples {M at n H ^ C?) 
with a well defined limit given by the triple (C°°(S^), Ti, C^). The number of degrees of freedom 
in the function algebra of is N"^ and the noncommutativity parameter is ^ = , . 



3 Theoretical predictions 



3.1 Gauge action 

It has been shown in [26-29] that the differential calculus on the fuzzy sphere is three 
dimensional and as a consequence, treating gauge fields as one-forms, a generic gauge field, A 
has 3 components. Each component A^, a = 1,2, 3, is an element of Mat]\f and the U{1) gauge 
symmetry of the commutative sphere will become a U {N) gauge symmetry on the fuzzy sphere 
with gauge transformations implemented as Aa — >UAaU^ -\- [/[!/„, C/^] where U e U{N). In 
this approach to gauge fields on the fuzzy sphere, S^, it is difficult to split the vector field A in 
a gauge-covariant fashion into a tangent gauge field and a normal scalar field. However, we can 
write a gauge-covariant expression for the normal scalar field as $ = ^{xaAa+AaXa+^=)- In the 
commutative limit, iV — cxo, we have A^ —>■ Aa and = UaAa and the splitting into gauge 

field and scalar field becomes trivial being implemented by simply writing Aa = Ua^ + aa, with 
nada = 0, where n is the unit vector on , (p = n ■ A is the normal gauge-invariant component 
of A and a is the tangent gauge field. 

For this formulation of gauge field theory on the fuzzy sphere the most general action (up 
to quartic power in Aa) on is then 



'S'Ar[^] — j^TrF"^^ — ^^^abcTr ^FabA^ — ^[Aa,Ah\A, 



(3.1) 



In above Trl = N and Fab = i[La, Ab\ —i[Lb, Aa\ + eabcAc + i[Aa, Ab\ is the covariant curvature, 
and Sn[^] = 0. The limit rn? — >oo gives a large mass to the scalar component and effectively 
projects it out of the spectrum of small fiuctuations. 

The associated continuum action 5*00 is then at most quadratic in the field and as a con- 
sequence the theory is largely trivial, consisting of a gauge field and a scalar field that have a 
mixing in their joint propagator. Indeed we can show that 



1 



An 



{fab? - ^eabJabricf - 2(£a¥')' + 4(1 + 2m^)^^ + Apip 



(3.2) 



where fab is the curvature of the tangent field and fab = ii^aO-b ~ ii^bO-a + ^abcO-c- As one can 
immediately see this theory consists of a 2-component gauge field that mixes with a scalar 
field ip, i.e. the propagator mixes the two fields. In the following we will primarily be interested 
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in the case with p = 0. We see also that the presence of the scalar field means that the geometry 
is completely specified, in that all the ingredients of the spectral triple are supplied by this field. 
In contrast a two dimensional gauge theory on its own would not be sufficient to specify the 
geometry. 



3.2 Matrix model 



We introduce a 



/? where (3 can be interpreted as an inverse temperature and we 



can rewrite the above gauge action (13. ip (shifted by constants and dropping the subscript A^) 
in terms of Da = La + Aa as follows: 



N 



2g'^N 



eabcTr 



1 i 

-FabA,- -[Aa,Ah]A, 

2 



~ 4 

a C2 



^Tr[Da, A]' + jeatcTrDaDbD, 



(3.3) 



V[D] 



2m2 
N 



m 
2^ 



■Tr{Dlf + 



/VC2 

2g^ 



m 



(3.4) 



The complete action functional is then: 



S[D] = + V[D] = Sn[A] - ^ 



-4 2 

2 



a m C2 I 



(3.5) 



and the constants are chosen so that ^[O] =0. The action takes the rather simple form 



S[D] = f Tr [ - \[Da, Dbf + feabcDaDbD, + g(I^D' - pD 



(3.6) 



where /x 



m 



. It is invariant under unitary transformations U (N) and global rotations 



5*0(3). Extrema of the model are given by the reducible representations of SU{2) and commut- 
ing matrices. For sufficiently small p and with C2 = (A^^ — l)/4 the classical absolute minima of 
the model is given by the irreducible representation of SU (2) of dimension A^. Small fluctua- 
tions around this background can then be seen to have the geometrical content of a Yang- Mills 
and scalar multiplet on a background fuzzy sphere as described in the previous section. As 
we will see, for small enough coupling or low temperature, these configurations also give the 
ground state of the fluctuating system. For very small negative values of the parameter /i, and 
for m = there is a local minimum at Da = and a global minimum at Da ~ La, separated by 
a barrier. As /i is made more negative the difference in energy (or Euclidean action) between 
the two extrema becomes less and eventually for /i = — | the two minima become degenerate. 



one occurring at Da = while the other occurs at Da 
barrier whose maximum occurs at Da = ^La. 



\La and they are separated by a 



In this special case (i.e. m — and // = — | or equivalently p — action takes the 

form 



a'^/i 1 \ 

S[D] = - [^Tr[Da, D,] + -6„,,D J . (3.7) 

and wc sec that the configurations Da = and Da = "^La both give zero action, however there 
is a unique configuration with Da = while there is an entire SU (N) manifold of configurations 
Da = '^ULaW which are equivalent. The classical model has a first order transition at m = 
for /X = — I and the classical ground state switches from Da — (pLa, where = {1 + ^/1 + 4//)/2, 
for // > — 2/9 to = for // < —2/9. The quantity (f) is therefore a useful order parameter 
for this transition. However, one would expect fluctuations to have a signiflcant effect on this 
classical picture. 

In fact both theoretical and numerical studies show that in the fluctuating theory, the fuzzy 
sphere phase only exists for p > —2/9 [30]. 

3.3 Quantization and Observables for small 

The quantum version of the model is taken to be that obtained by functional integration 
with respect to the gauge field. This amounts to integration over the three Hermitian matrices 
Da with Dyson measure and the partition function Z is given by 

Z = J dDae-'^^^ = e-^^^^s"' J c^X„e-^^[^] , X„ = a = aVN 

S[X] = -\Tr[Xa, + ^-^ea,cTrXaX,X^ + "^TriXlf + [^0= - mW^ Tr{Xl). 

(3.8) 

The latter form of the expression (in terms of Xa) allows us to take a = and we see that this 
limit is equivalent to removing all but the leading commutator squared term. Also the quartic 
term proportional to m? survives. However if m? and p are scaled appropriately with a only 
the Chern-Simons term eabcTrXaX^Xc is removed. 

The set of gauge equivalent configurations is parameterized by the SU (N) group manifold 
which is compact, so there is no need to gauge fix and the functional integral is well defined, 
being an ordinary integral over R^^^ . However, the volume of the gauge group diverges in the 
limit N ^ oo and to make contact with the commutative formulation it is convenient to gauge 
fix in the standard way. 

In the background field gauge formulation we separate the field as Xa = ctDa + Qa- The 
action is invariant under Da — >Da, Qa — ^UQaU"^ + U[Da, U^. 

Following the standard Faddeev-Popov procedure [31] and taking the background field con- 
figuration to be Da — (f)La one finds, keeping fixed as the — > oo limit is taken, that 
F = — In Z is given by 



^ = |logci^ + ^ 



4 3 ^ 4 ^2 



+ log q;0. 



(3.9) 
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m 



with /i 

The most notable feature of this expression is that the entire fluctuation contribution is 
summarized in the logarithmic term log(a;0). From (13.91) we see that the effective potential for 
the order parameter (j) is 



2C2 



\- m u — 

4 3 4^2 



+ log< 



(3.10) 



The effective potential is not bounded below at = due to the ln0 term. However, our 
analysis assumes the existence of a fuzzy sphere ground state and so the effective potential can 
only be trusted in this phase. It has a local minimum for (p positive and sufficiently large a and 
in this regime the fuzzy sphere configuration exists, for lower values of a our numerical study 
indicates that the model is indeed in a different phase. 
Setting (j)^ = (or equivalently (j^^^-^r^ = 0) gives 



+ 1 = 0, 



(3.11) 



the solution of which specifies 0. We define the average of the action, which will be one of the 
principal observables of our numerical study, as 



S=<S> /N^ 



a' 



Then a direct computation and use of eq. (13. lip yields 



24 



We can also compute the expected radius of the sphere 



1 

R 



< TrPj > 



2 dF 
dfi 



This can also be calculated directly in perturbation theory as 



(3.12) 



(3.13) 



(3.14) 



1 
R 



+ 



;Tr3TrN^{- 



1 



Hence, we conclude that the expected inverse radius of the fuzzy sphere is given by 



Also 
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(3.15) 



(3.16) 



(3.17) 



which yields 



<Tr{Dlf> ^ ^ 
2Af3c2 8 • 



(3.18) 



Scahng {Da)ij to {l—e){Da)ij in both the action and measure amounts to a simple coordinate 
transformation which leaves the partition function invariant. However it also leads to the 
nontrivial identity 



2m^ 

< >= 3N'' , Km = Tr{ - [D„ D,]'' + 2ieabcDaD,D, - 2m^Dl + (Dl)'' ) . (3.19) 



Using this identity we can express 



C2 



^ 4 ^ 6A^3 
or equivalently in the form 



< ieabcTrDaDbDc > 



a 



2m 



< TrDl > . 



(3.20) 



5 = 1 + 



a 



<Tr[Da,Dbf > 



a 



W < TrDt > 



m' < Tr{Dlf > 



We define the Yang-Mills and Chern-Simons actions by 



4YM = -<^:iM£> , and 3CS = <i^abcTrD.D,D.> _ 

2Nc2 ' Nc2 



Then using the above results we find 



4YM 



+ — and 3CS 

a'* 



(3.21) 



(3.22) 



(3.23) 



In the above we have extensively used the fact that must satisfy (13.111) and taken p = i.e. 

Another significant observable for our numerical study is the specific heat Cy defined as 

_<S^>-<S>^ 



ATS 



< S > d i < S > 

— a 



iV2 



A direct calculation yields 



- + 

4 ^ 32 



(3.24) 



(3.25) 



Finally one can recover perturbation theory in the coupling (7^ = 1/a^ by expanding in 
In particular the one-loop predictions are obtained by using the solution of ( 13. lip expanded to 
first order which is a — *^cxo given by 



(3.26) 



In the next section we will look at the solution of (13.111) in more detail and the consequences 
for the transition. 
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3.4 Phase Transitions 



Let us now examine the predictions for the quantum transition as determined by F given 
in (13.91) or the effective potential (13.101) which we repeat here for convenience 



2c2 



+ log. 



(3.27) 



But first let us review the classical case. The only difference between the full quantum potential 
(I3.10p and the corresponding classical potential is the quantum induced logarithm of 0, which 
as we will see plays a crucial role. The extrema of the classical potential occur at 



;i + ,n')4. = { 0. 1±4±5 j . (3,28) 



where t = /i(l + m^). The first and last expressions are local minima and the middle one is 
the maximum of the barrier between them. For fi positive the global minimum is the third 
expression, i.e. the largest value of 0. When written in terms of we see the potential takes the 
same form as that for m = and we can read off that if /i is sent negative then this minimum 
becomes degenerate with that at = at t = /i(l H-m^) = — | and the maximum height of the 
barrier is given by 324(1+^^)3 • When /i = the minimum is clearly 0=1, i.e Da = La for all 
m and this is separated from the local minimum at = by a barrier. The classical transition 
therefore has the same character as that of the m = model and the transition is 1st order and 
occurs only when p is tuned to a critical value. 

Let us now consider the effect of the fluctuation induced logarithm of 0. The potential is 
plotted in figure [T] for different values of a and fi = m? = 20. The condition V^ff{(j)) = gives 
us extrema of the model. For large enough a (or low enough temperature) and large enough 
m and /x it admits four real solutions two for positive and two for negative 0. The largest of 
the positive solutions can be identified with the least free energy and therefore the ground 
state of the system in this phase of the theory. It will determine the actual radius of the sphere. 
The second positive solution is the local maximum (figured]) of Vcs{4>) and will determine the 
height of the barrier in the effective potential. As the coupling is decreased (or the temperature 
increased) these two solutions merge and the barrier disappears. This is the critical point of 
the model and it has no classical counterpart since, in the classical case, the barrier between 
the two minima never disappears. For smaller couplings than this critical coupling a^, the fuzzy 
sphere solution Da = (pLa no longer exists and the effective potential cannot be relied on. This 
is in accord with our numerical simulations which indicate that as the matrix size is increased 
the radius as defined in (13.141) appears to go to zero. 

The condition when the barrier disappears is 

~ 4 1 
Kff = y [302 - 20 + 3m20' - /x] - = 0. (3.29) 
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Solving both fl3?TT|) and iK29h yields 4(1 + m^)^^ - 30* - 2/i = and therefore the critical 
values 



1(1 + 



1 + 



-2/x) 



(3.30) 



(3.31) 



where as defined earlier 
yields 



0/(1 + m^) and t = /i(l + m^). Setting both = and = 



3 

4 ' 



alone leads to no significant simplification and we still have 



(3.32) 
with 



while setting m 
t = fi. 

If we take fi negative as in the discussion of the last section we see that g^, goes to zero at 
t = —1/4 and the critical coupling a^, is sent to infinity and therefore for t < the model has 
no fuzzy sphere phase. This case arises when 

1 



4(1 



(3.33) 



However, in the region — ^ < yU < — | the action (13.61) is completely positive. It is therefore not 
sufficient to consider only the configuration Da = (pLa, but rather all SU (2) representations 
must be considered. Furthermore for large a the ground state will be dominated by those 
representations with the smallest Casimir. This means that there is no fuzzy sphere solution 
for /i < — |. A result that we also observe in simulations and in agreement with the result 
of [30]. We therefore see that the classical transition described above is significantly affected 
by fluctuations and in particular the fuzzy sphere phase dissapears when p is increased to the 



special value of 



The other limit of interest is the limit p 



m 



>oo. In this case 



1 

V2 



at 



(3.34) 



This means that the phase transition is located at a smaller value of the coupling constant a as 
m is increased. In other words the region where the fuzzy sphere is stable is extended to lower 
values of the coupling or higher temperatures. These results agree nicely with numerical data. 

As we cross the critical value of a, a rather exotic phase transition occurs where the geom- 
etry disappears as the temperature is increased. The fuzzy sphere phase has the background 
geometry of a two dimensional spherical non-commutative manifold which macroscopically be- 
comes a standard commutative sphere for oo. The fluctuations are then of a U{1) gauge 
theory which mixes with a scalar fleld on this background. In the high temperature phase, 
which we call a matrix phase, the order parameter is not well deflned and the fluctuations are 
around diagonal matrices so the model is a pure matrix one corresponding to a zero dimensional 
Yang- Mills theory in the large A^ limit. The fuzzy sphere phase occurs for a > a^, while the 
matrix phase occurs for a < a*. 
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3.5 Predictions from the effective potential 



Since = is not a solution of equation fl3.1ip the extremal equation can be rearranged 
and when expressed in terms of (j) takes the form 



pa) = 04 _ 73 _ ^72 ^ 



0. 



(3.35) 



where = a^/(l + m"^)^. By the substitution = (1 — Ax)(f) + x, P{(f)) can be brought to the 
form of a 



- 0^ - + ^ = 0. 



(3.36) 



The new parameters A and f3 are given by 

X 



A 



r 

;i _4x)3L^+ 3 J 



2(l + m2)3 1 



«4 



m^(l + m^l 



/34 (1 - 4x) 
and the shift x must take one of the two values 

1 



X + 



5m (1 + m 



(3.37) 
(3.38) 



X^X± = -[l±^l + j). 



(3.39) 



We choose x = x_ since x_(0) = and so this case allows us to easily recover the case with 
H = = and for — >0 we get 



x_ = -— + 0(m*) , A = O(m^) 
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— (1 - -m^ + 0(m^)) 



3 ■ ' ' ' ' /34 a4^" 3 

The two positive solutions of (13.351) . (f)±, for general values of and fi are then given by 



(3.40) 



1 + (1 - 4x)Vl + 5±(l - 4a;) J2 - 5 + 



(1 + 8A) 



5 = 4Vr3 



A 



(3* 2 



(3.41) 



V3/34 3 

We rewrite the above solution ( ]3.4ip as follows, (collecting definitions here for completeness) 



t = jj{m'^ + 1) , 0-^ 



:i + m2)3 



, 8t aH^ 

13 



(§. _|_ allies 
2 V3 ~'~ 9 ^ 



(3.42) 



Then we can show that 

(1 -4a;)^(l + 8A) = 1 + 4t , W 

and 



1 V = l-^ 



Substituting one finds the relatively simple form 



h± = (1 + m2)0± - ^ , 2 4 



3 2-^ 2 3 4^i + f +d 



For completeness the remaining two solutions of the quartic P{4>) = are given by 



ij = (l+m2)0- = 1- 1 /i + f + rf±i 



2A/ 2 



+ 



l+4t 



(3.43) 



(3.44) 



(3.45) 



(3.46) 



At the critical point a = p becomes zero and the two solutions 0+ and are equal. 

In the fuzzy sphere phase the ground state of the system is given by 0+ and the barrier 



maximum is at 0_. The minimum 0+ together with the powers 



and (jA are plotted in 



figures [2] and [3] for w? = and 200. 

In the specific heat we will need the derivative of the minimum with respect to a. This is 
given by 



da 



1 dd 
4 da 



1 + At 



1 



l+4t 



' + 'i+d) 



+ 



4y|Tf^^4 3 
1 



1 I 2t 



d- 



(3.47) 



dd _ 4(1 + m2)^a 
da 3a^ 

4(1 + m2)^at 
355 



1 I 8t ^ (l + f)g 4(| + i^ 
3 ^ 



1,8^ a + f)g , 4(i + ^ 

3 ^ 



1 

(g- ^)i 



(3.48) 



3.6 Critical behaviour 

For fi = = we have x = x_ = and A = and the solution 



1+ VI+5+ ,2-5+ 



_ 4 
as 



i + >/i-§)Vii 



1 - 



at 



(3.49) 
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Expanding near the critical point we have S = 3 — ye + 0(e^) and we obtain thus the expression 



6e - f + 0(e 



(3.50) 



Substituting into (I3.13P near the critical point we obtain the expression for the scaled average 
action 



S 



12 



r^V" - "* r^{a - a^) + 0{{a - a^y^) , 



3H2H 



and the specific heat is then given by 



29 
36 



+ u , - + 0{{a - a*)^). 



This gives a divergent specific heat with critical exponent 



a 



2 ■ 



(3.51) 



(3.52) 



(3.53) 



If instead we consider fi = but m 7^ so that t = then the above expressions remain 
essentially the same provided we substitute a for a. The critical behaviour remains the same 
for this case with the critical inverse temperature 



(3. = at = qr{l+my 



(3.54) 



so that increasing with /i = sends the critical temperature lower and so the region of 
stability of the fuzzy sphere solution is reduced. 



3.7 The generic case 

At the critical point the coupling a takes the value a^, and the order parameter takes 
the value 0=,,. At this critical point the two solutions and 0_ merge. This is more easily 
determined by requiring the additional equation 



2(l+m2)3 dVeff 



g(0) = 30^ - 203 _ t02 - — = 0. 



Putting the two equations -P(0) =0 and Q{4>) = together we obtain 

-a -jn4>) 



P(0)+Q(0)=40^- 



2t0" 



(3.55) 



(3.56) 



In other words 



dP{i) 



at the critical point. Expanding around the critical 

= 0* + a. 
15 



(3.57) 



and using P{(t)*) = P'(0*) = we obtain 



1W + 



at 



0. 



(3.58) 



For small a, treating cr^ and perturbatively we get 



(40, - 1) f]__]_ 

A A \ 



30, + 4tVa* (30, + 4t)2Va4 

1 , (160, - 1) 



xl ^30* + 4t 



a — a, — 



a,5(30, + 4t)^ 



(a - a,) + 



Hence 



+ \ ~ + ••• 



or equivalently to leading order we have 



\ 5 



i=Va - a, + ... 



(3.59) 



(3.60) 



(3.61) 



The average action S near the critical point can be computed using this expression of in 
equation (I3.13p . The result is that 



4 0* (0* + /« - «c 

o — o, — d 



Where 



i5. 



^30, + 4t V a, 

_ 3 _ (0, + 3t) 
" 4 3(0, + 2t) 



+ ... 



(3.62) 



(3.63) 



which interpolates between 5, = ^ for t = and iS, = 1 for large t. 
In order to compute the specific heat we need the derivative 



da 



^1 a/30, + 4/i a/o - a. 



+ ... 



(3.64) 



The divergent term in the specific heat is still given by a square root singularity. From (13.251) 

we get 



+ ... 



where the background constant contribution to the specific heat is given by 

5^3 (3 + 4t)0, + 2t 
" 4 8(30 + 4t)2 
16 



(3.65) 



(3.66) 



If we set ;U = in (13.611) and (13.651) we recover (I3.50p and (I3.52p . Let us recall that the critical 
value a* can be given by the formula 



at 



Then (13.651) can be put in the form 



Cv — C., + 



1 



8\/l+ 



+ ... 



(3.67) 



(3.68) 



The prediction here is that the critical exponent of the specific heat for this model is given 
precisely by 



1 

"=2 



(3.69) 



Specializing to the case /i = we see the coefficient of the singularity for any small (i.e 
the amplitude) is 



c(m^) 



1 + 



54 03 
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If we extrapolate these results to large where we know that 
we get 



1/^2 and at 



Cv 



(3.70) 

> 8/m2 
(3.71) 



The coefficient of the singularity and the critical value a^, become very small and vanish when 



m 



oo. For comparative purposes we can compute the ratio 



c(m J 

7(0)" 



-i-Y- 



(3.72) 



For = 200 we get the ratio c(200)/c(0) = 0.57 which is not yet very small. As we will see 
below our data (see figure [T2!) in the critical region for large m is not precise enough to confirm 
or rule out the presence of a singularity. 



4 Numerical results 

Let us now turn to the numerical simulations. A fully nonperturbative study of this model 
is done in [7, 15]. For the model with = see also [5, 17]. In Monte Carlo simulations we 
use the Metropolis algorithm and the action (13.51) with in the range A^ = 8 to 104 and in 
the range = to 1000. The errors were estimated using the binning-jackknife method. We 
measure the radius of the sphere R (the order parameter) defined by NC2/R =< TrD^ >, the 
average value of the action < S > and the specific heat = '^'•'^"^f^-* ^ as functions of a for 
different values of A^ and m^. We also measure the eigenvalue distributions of several operators. 
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Figure 1: The effective potential for = 20. The minimum corresponds with the fuzzy sphere 
solution, as we lower the coupling constant the minimum disappears and the fuzzy sphere 
collapses. 







8 10 



Figure 2: The functions 0, 0^, —0^ and (j)'^ for = 0. corresponds with 0+; the minimum 
of the effective potential (eq l3.27p . 
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Figure 3: The functions 



—d)^ and d)'^ for m? = 200. 



4.1 The theory with = 

For = we observe that the expectation values S, YM, CS (defined above) are all 
discontinuous at as = 2.1 ±0.1 (figureHj). This is where the transition occurs. Indeed this agrees 
with the theoretical value a^, = 2.087. The theory also predicts the behaviour of these actions 
in the fuzzy sphere phase. There clearly exists a latent heat and hence we are dealing with a 
1st order transition which terminates at some value of m^. The radius is also discontinuous 
at the critical point (figure [5]) whereas the specific heat is discontinuous and divergent (figure 
|6]). Near the critical point we compute a divergent specific heat with critical exponent a = 1/2 
(equation (13.521) ). The fit for data fixing a = 0.5 gives the critical value ac = 2.125 ± 0.007 
again with good agreement with the theory. From the matrix side the specific heat seems to be 
a constant equal to 0.75 and therefore the critical exponent is zero. The inverse radius in the 
matrix phase goes through a minimum and then rise quickly and sharply to infinity. 



fuzzy sphere («>«=,< ) 


matrix phase (a < a*) 


R = 1 


R = 


a = l 


Cy = 0.75 
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Figure 4: The observables CS, YM and ^ for = function of the couphng constant 
for different matrix sizes N. The sohd hne corresponds to the theoretical prediction using the 
local minimum (13.451) of the effective potential. 
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Figure 5: The inverse radius for = 0.25 against a for N = 10, 16, 24. The solid line 
corresponds to 0^, given by eg. (13.451) . 
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Figure 6: The specific heat for = as a function of the couphng constant for = 
16, 24, 32,48. The curve corresponds with the theoretical prediction given by eq. fl3.25p for 
= 0. 

4.2 The action, radius and specific heat for ^ 

For small values of we determine the critical value ag as the point of discontinuity in 
5,YM, CS, Tr{Diy and the radius TrD^. This is where the divergence in C„ occurs. For 
example for = 0.75 the action looks continuous but its parts are all discontinuous with a 
jump. The radius is also discontinuous with a jump. The specific heat is still divergent in this 
case (figure [7]). This is still 1st order. 

For = 55 the action and its parts become continuous. We find in particular that the 
Chern-Simons and the radius are becoming continuous around = 40 — 50. These are the 
two operators which are associated with the geometry. The specific heat seems now to be 
continuous (figure [H]). This looks like a 3rd order transition. 

4.3 The Umit — > oo and specific heat 

In this case we measure the critical value ag as follows. We observe that different actions 
< S > which correspond to different values of (for some fixed value of m^) intersect at some 
value of the coupling constant a which we define as (figure [H])- This is the critical point. For 
example we find for = 200 the result = 0.4 ±0.1. The theoretical value is a^, = 0.44. For 
large the theoretical critical value a* is given by equation fl3.34p . The measured value cts 
tends to be smaller than this predicted value. 

The quantities S, YM, CS, {Tr{Dl)^) and the radius are all continuous across the tran- 
sition point in this regime (figure [TUl) . We observe that near the critical point the numerical 
results approach the theoretical curves as we increase N. We also checked the Ward identities 
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Different observables for = 0.75 plotted against a for different matrix sizes. Again 
lines represent the theoretical predictions using the local minimum of the effective 
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Figure 9: The observables (left) and < S > (right) for = 200 plotted as functions of a 
for N = 16, 24, 32, 64, 104. The value of a at which the curves < > for different values of N 
cross is defined as a,. 



( 13A9D . ( 13:2011 and (^23) (figure [III) . 

The specific heat in the fuzzy sphere phase is constant equal to 1, it starts to decrease at 
<5max, goes through a minimum at amin and then goes up again to the value 0.75 when a — > 
(figure [12]). The values <5max,min decrease with while the minimum value of increases. 
Extrapolating the amax and amin to = oo (figure [131) we obtain our estimate for the critical 
coupling which agree with within errors. The matrix-to-S"^ phase transition looks then 
3rd order. However it could be that for large N the specific heat becomes discontinuous at the 
critical point with a jump, i.e the transition is discontinuous with 2nd order fluctuations. 

Thus it seems that the specific heat in the regime of large values of is such that i) 
fimax.min approach ctg in the limit — > oo and ii) that the specific heat becomes constant in 
the matrix phase and equal to = 0.75. There remains the question of whether or not the 
specific heat has critical fluctuations at the critical point for large m^. 

The theory still predicts a transition (equation (I3.65P ) with critical fluctuations and a di- 
vergence in the specific heat with a critical exponent a = 1/2 and with a very small amplitude 
(the coefficient of the square root singularity). There is possibly some evidence for this even for 

= 100 but none for = 200. To resolve this question we need to go very near the critical 
point and simulate with bigger A^. 

4.4 The phase diagram 

Our phase diagram in terms of the parameters a and which we have studied is given in 
figure [HI We have identified two different phases of the matrix model 03.51) . In the geometrical 
or fuzzy sphere phase we have a U{1) gauge theory on Sf^; the geometry of the sphere and the 
structure of the U{1) gauge group are stable under quantum fluctuations so the theory in the 
continuum limit is an ordinary U{1) on the sphere. In the matrix phase the fuzzy sphere vacuum 
collapses under quantum fluctuations and there is no underlying sphere in the continuum large 
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Figure 10: Other observables for = 200 as functions of the couphng constant for = 
16, 24, 32, 64. The sohd hues correspond to the theoretical predictions. We observe that the 
data tend to approach the theoretical prediction as is increased. 
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Figure 11: The Ward identity (13.191) for = 200 as a function of the coupling constant a for 
N = 10,16,24,32. 
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Figure 12: The specific heat for = 100 and 200 for different matrix sizes. In the region of 
the fuzzy sphere the value of the specific heat is 1. In the region of the matrix phase the value 
of the specific heat tends to take the constant value 0.75 as the size of matrix increases. 
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limit. In this phase the model should be described by a pure matrix model without any 
background spacetime geometry. The transition in the Ehrenfest classification would be labeled 
a first order transition, but this classification is not very helpful. The transition described here 
is a very exotic one with both a latent heat and a divergent specific heat. We know of no 
other example of such a transition. As we follow this line of transitions the entropy jump 
or equivalently the latent heat becomes zero at around ~ 40 and remains zero for larger 
m^. Our theoretical analysis indicates that there is still a divergent specific heat, however our 
numerical simulations are not fine enough to determine whether this is so or not. 

For large the transition (we expect) still has a divergent specific heat as the transition 
is approached from the fuzzy sphere side. This is the conclusion of our theoretical analysis and 
our numerical results are consistent with this. But our numerical results are not conclusive. It 
may also be that the transition is even in Ehernfest's classification a 3rd order, where there 
is a jump in the specific heat with no divergence. We could not determine the nature of the 
transition in the large m regime with any confidence from the numerical data. In all cases 
the fuzzy sphere-to-matrix theory transitions are from a one-cut phase (the matrix phase) to a 
point or discrete spectrum in the geometrical (fuzzy sphere) phase. 

The specific heat in the fuzzy sphere phase takes the value 1 where the gauge field can only 
contribute the amount 1/2. This can be understood as follows. The high temperature limit 
of the specific heat for any matrix model is governed by the largest term and must go like 
Ntotai I degree where Ntotai is the total number of degrees of freedom and degree is the degree 
of the polynomial. This gives the limiting high temperature limit of the specific heat be | for 
all values of the parameters. In the simple model with m = and = this value is achieved 
from the transition point onwards. 

For the full model in the large m? regime the effect of the potential V should be dominant. 
This can be seen by remarking that in the strong-coupling limit m? — >oo, a — >0 keeping fixed 
a^m^ the action S reduces to V . Note that if we consider V alone with a measure given by 
(with $' = \J D1) instead of j[dDo\ then we will get the usual quartic potential dynamics with 
a well known 3rd order transition. Here when we consider the model given by the potential 
V with the measure j[dDa\ we obtain the specific heat given in figure [151 In the region of 
parameters corresponding to the matrix phase the specific heat shows in this case a structure 
similar to that of the full model S. However in the region of parameters corresponding to the 
fuzzy sphere phase the specific heat is given now by C„ = 1/2. Thus the field $ contributes only 
the amount 1/2 to C„. Indeed from the eigenvalue distribution of the operator $ computed in 
the fuzzy sphere phase it is shown explicitly that $ has Gaussian fiuctuations. The behaviour 
of Cy in the full model S is thus a non-trivial mixture of the behaviours in Sq (first order 
transition) and V (3rd order transition) considered separately. 

Given that the behaviour of the full model can be described as a non-trivial mixture of the 
model 5*0 and the potential we expect that the effect of adding the potential to the model is 
to shift the transition temperature and provide a non-trivial background specific heat. The 
divergence of the specific heat arises from the interplay of the two terms in 5*0, i.e. between 
the the Chern-Simons and Yang-Mills terms. Given that this competition leads to a divergence 
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Figure 14: The phase diagram shows the curve separating the geometrical and matrix phases 



of the model (13. 6p with /i 



m 



The critical curve is given by eg. (13.311) and describes (at 



least for small values of m^) a line of exotic transition with a jump in the entropy yet with 
divergent critical fluctuations and a divergent specific heat with critical exponent a = 0.5, when 
approached from the fuzzy sphere side. The points which lay off the critical line in the middle 
mass region might suggest the existence of a multicritical point. 



of the specific heat it should eventually emerge from the background sufficiently close to the 
transition. 

4.5 The eigenvalue distributions for large 

4.5.1 The low temperature phase (fuzzy sphere) 

Numerically we can check that the normal scalar field and the tangent gauge field decouple 
from each other in the "fuzzy sphere phase" in the limit — >oo by computing the eigenvalues 
of the operators and D^. See figure [IHl The fit for the distribution of eigenvalues of — C2 
is given by the Wigner semi-circle law 



p(^) = zT^y^ls-^^- (4-1) 



This distribution is consistent with the effective Gaussian potential 



Asphe _ 2 „ \2 2 4:C2N 



= zrTriDi - c,Y , ai^ = (4.2) 



We find numerically 
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Figure 15: The specific heat for the pure potential with = 200 as function of the couphng 
constant computed with respect to the Boltzmann weight e~^. 
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The parameter a is the theoretical prediction given by = 4c2N/m'^a'^ which goes hke A^'^. 
The effective parameter a^g is found to behave as 

al^ = (3.7645 ± 0.6557) x 10"^ x iv3-3i70±o.o492_ ^4 3) 

This means that the renormahzed value aes of the gauge coupling constant is slowly decreasing 
with A^. Equivalently the parameter Oeff yields a small correction to the classical potential V 
which is linear in $. Indeed we can show that 



Z = / dDaC "off 



Co = -2Ncl(--—) -^NHog^. (4.4) 

This shows explicitly that having 7^ a means that there is an extra linear term in added 
to the classical potential. 

In the fuzzy sphere phase the field configurations are thus given by (or are close to) 
representations of SU{2) of spin s = as we can clearly see on figure [161 Indeed the eigen- 
values of for a = 5 and m? = 200 are found to lie within the range —^^y^, ■ ■ ■ , 0, ■ ■ ■ , ^^-^ 
as expected for iV = 24 and = 32. The eigenvalues of the commutator —i[Di, D2] are also 
found to lie in the range ■ ■ ■ , 0, ■ ■ ■ , ^^^^ (figure fTTll . 
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Figure 16: The eigenvalue distribution of — C2 (left) for = 24,32,48 and (right) for 
= 24, = 32 and = 48 in the fuzzy sphere phase with = 200 and a = 5. The fit for 
Da — C2 corresponds with the Wigner semi-circle law (14. ip . 
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4.5.2 The high temperature phase (matrix phase) 



The distribution of the eigenvalues of suffers a distortion as soon as the system un- 
dergoes the matrix-to-sphere transition and deviations from the Wigner semi-circle eigenvalue 
distribution (14.11) become large as we lower the coupling constant a. See figures [T8] and [T9l 

In this phase the distribution for D3 is symmetric around zero and the fit is given by the 
one-cut solution 



p{x) 



16 



c(c + 46) 27r 



(x^ + b) vc — a?. 



(4.5) 



By rotational invariance the eigenvalues of the other two matrices Di and D2 must be similarly 
distributed. This means in particular that the model as a whole behaves in the "matrix phase" 
as a system of 3 decoupled 1-matrix models given by the effective potentials {i fixed) 



V^matr 
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'eff 



/2N _ 3c\ 



TrDf + -^TrD' 
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^eff 



c2 2cb 



4:C2N 
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eff 



We find numerically that 



N 


b 


c 


24 


1047.23 ±15.97 


633.655 ±0.7516 


32 


2392.16 ±65.09 


1074.07 ± 1.8680 


48 


10350.9 ± 236.70 


2185.85 ± 1.7130 



Again the theoretical prediction for goes like A^^ whereas the effective parameter ""^ 
to behave as 



is found 



(0.6147 ±0.3059) x iV' 



3.6580±0.1553 



(4.7) 



matr 
eff 



By going through the same argument which lead to equation (14. 4p we can show that 
is equivalent to the addition of an extra linear correction in $ (which depends on c and b or 
equivalently c and a^g) to the classical potential 



4c2 9 2 4 
—TrD^ H TrD^. 

a2 ' ^ a? ' 



(4. 



The result (14. 6 p accounts for the value = 0.75 of the specific heat. Indeed for effective 
potentials of the form (14.61) each matrix Di contributes the amount 0.25 . 

A final remark is to note that the eigenvalue distributions for — C2 in figure [16] and in 
figure dn] clearly depend on A^. It would be desirable to find the proper scaling of the parameters 
for which the distributions are A^-independent. This is also related with the predictions of the 
effective parameters and the corresponding effective potentials written in (14.21) and (14. 6p . 
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the fit ( 14. Sp which is the the eigenvalue distribution in the regime of one-cut solution. 
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4.6 Emergent geometry 

The sequence of eigenvalue distributions for a = 0.2,0.6, 1, 1.5,2 and 4 for = 200 and 
iV = 24 (figure [20]) show clearly that a geometrical phase is emerging as the coupling is increased 
or equivalently as the temperature is reduced. The geometry that emerges here is that of the 
fuzzy sphere. This geometry becomes the classical sphere as is sent to infinity. This is the 
geometry of the sphere emerging from a pure matrix model and is in our opinion a very simple 
demonstration of a novel concept and opens up the possibility of discussing emergent geometry 
in a dynamical and statistical mechanical setting. We see clearly that as a is increased from a 
value deep in the matrix phase to a value well into the fuzzy sphere phase that goes, from a 
random matrix with a continuous eigenvalue distribution centered around 0, to a matrix whose 
eigenvalues are sharply concentrated on the eigenvalues of the rotation generator L3. This is 
also true for the matrices Di and D2 which go over to Li and L2 respectively in the fuzzy 
sphere phase. This can be seen explicitly in figure [T7I since the commutator [Di,D2] is found 
to be well approximated by the matrix iL^. 

4.7 The model with only pure potential term 

The model is given by the potential 



The configurations which minimize this potential are matrices which satisfy = C2- We plot 
(figure [?T]) the eigenvalue distributions of the matrices and D"^ for = 24, = 200 and 
a = 0.4, 0.6, 1.4. The distributions of seem to behave in the same way as the distributions 
of D"^ computed in the full model for all values of a. The distributions of D3 are similar to the 
corresponding distributions of computed in the full model only in the region of the matrix 
phase where they are found to fit to the one-cut solution (14.51) . In the region of the fuzzy 
sphere phase the distributions of D3 tend to split into two cuts. However they never achieve 
this splitting completely due to the mixing terms TrDfDl, TrD\D\ and TrDlD^. These terms 
cannot be neglected for values of the coupling a where the model is in the fuzzy sphere phase. 

4.8 The C hern- Simons+ potential model 

In this final section we present the eigenvalue distributions for the model in which we set 
the Yang-Mills term to zero. The action reduces to 
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~ A 1 
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The equations of motion are 



teabc[Db, D,] - 2m^Da + —{Dl Da} 
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Figure 20: The model undergoes a transition from the matrix phase with no underlying geom- 
etry to a geometrical phase as the coupling constant a"^ = j; is increased (or the temperature 
lowered). The geometry which emerges in this model is the geometry of the sphere. 
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Figure 21: The eigenvalue distributions of and for the pure potential term. Deep inside 
the nonperturbative regime (small a or large temperature) the distribution for D3 is in the 
one-cut solution (14.51) whereas for large a (or low temperature) the eigenvalue distribution for 



is given by the Wigner semi-circle law. 



The possible solutions are the commuting matrices and reducible representations of SU{2). 

We measure the eigenvalue distributions of the matrix for different values of a. The 
distributions of D3 in the region of the matrix phase can be fit to the one-cut solution (14. 5p . 
In the region of the fuzzy sphere phase the distributions of split into two well separated 
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Figure 22: The eigenvalue distributions of for the Chern-Simons + Potential model for 
different values of the coupling a for N = 24 and = 200. In this case, the system undergoes 
a transition from the one-cut regime to the two-cut regime as the coupling (or temperature) is 
increased (lowered). 



cuts. See figure [221 This model exhibits therefore the behaviour of a typical quartic one- matrix 
model. The value of the coupling a where the transition from the one-cut phase to the two-cut 
phase happens coincides with the maximum of the specific heat. From the numerical results 
we can also see that in the regime of large a the preferred configurations are given by 

Da = 2X^®lN. (4.12) 

A is fixed by the requirement = C2, i.e 3A^ = C2 and aa are the Pauli matrices. Indeed the 
eigenvalue distributions of the matrices D^, DaDhDaDf,, DaDhDhDa and ieabcDaDbDc for large 
a are found (figure 1251) to be given by the Wigner semi circle laws 

2 

p(x) = ^^ya^ - (x - xo)^. (4-13) 
a 71 

The centers Xq are given by the theoretical values 

Dl = 3XHn 

DaDbDaDb = -3XHn 

DaD^DbDa = 9A^1^ 
leabcDaD^D, = -eA^l^. (4.14) 

For instance in figure [231 we can see from the eigenvalue distribution of D3 with = 24, a = 4 
and = 200 that the preferred configurations are given by (14.121) with A = ±(7.16 ± 0.25) 
whereas the predicted value is A = ±6.9567. We also remark that the eigenvalue distributions 
of Dl and the other operators have the same structure as in the full model in its matrix phase. 
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Figure 23: The eigenvalue distributions of D^, D^, DaDhDaDf,, DaDhDbDa, ieabcDaDbDc for 
the Chern-Simons+potential model. The solid lines correspond to the Wigner semicircle laws 
with centers given by eqs.f l4.14]) . 



5 Conclusion and outlook 



We have studied the simple three matrix model with Euclidean action functional (13.60 for 
general values of its parameters (3 = and m but focused on a small range of the possible 
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values of the parameter ji. 

We find the model to have two clearly distinct phases. In the high temperature regime (i.e. 
small a) the model has a disordered phase. In this phase the eigenvalues distribution of an 
individual matrix is well approximated by the one-cut distribution of a hermitian matrix model 
with quartic potential. We call this phase the matrix phase of the model. 

At low temperature the model has an ordered phase. The order is unusual in that it 
describes the condensation of a background geometry as a collective order of the matrix degrees 
of freedom. We call this the geometrical or fuzzy sphere phase, since for finite matrix size the 
ground state is described by a fuzzy sphere which is a quantized version of the classical sphere. 
In the large matrix size limit the sphere becomes classical but at a microscopic level the geometry 
always has a noncommutative character as can be seen from the spectrum of the "coordinate 
functions" which are proportional to the su{2) generator L3 in the irreducible representation 
of dimension given by the matrix size. 

In the geometrical phase small fluctuations arc those of a U{1) gauge field and a neutral 
scalar field fiuctuating on a round two sphere. The two fields have non-trivial mixing at the 
quadratic level but are otherwise not interacting. In this phase, with fi = rn^, the parameter rri^ 
parameterizes the mass of the scalar fluctuation, otherwise for /i — — p/N, the parameter p 
provides a constant external current for the scalar fleld. 

For m = fi = the transition between the two phases is found to be discontinuous. There 
is a jump in the internal energy (expectation of the Euclidean action) and from a theoretical 
analysis (valid in the fuzzy sphere phase) we find the entropy drops by 1/9 per degree of 
freedom [15] as one crosses from the high temperature matrix phase to the geometrical one. 
Our theoretical results also predict that: For all m and /x > — |, the model has divergent critical 
fluctuations in the speciflc heat characterised by the critical exponent a — 1/2. The critical 
regime narrows as the critical temperature decreases and the transition temperature is sent to 
zero at // = — I and so there is no geometrical phase beyond this point. 

Our numerical simulations are in excellent agreement with these theoretical predictions and 
we flnd the critical fluctuations arc only present in the fuzzy sphere phase so that the transition 
has an asymmetrical character. We know of no other physical setting that exhibits transitions of 
the type presented here and further numerical and theoretical study arc needed. However, the 
thermodynamic properties of the transition are similar to those found in the 6— vertex model 
and the dimer model [33,34]. 

We focus simulations on p, — m? and observe that the discontinuity in the entropy (or 
the latent heat) decreases as is increased with the jump vanishing for sufficiently large m. 
We have not been able to determine with precision where the transition becomes continuous, 
however the jump in the entropy becomes too small to measure beyond ~ 40. Also, the 
predicted critical fluctuations are not seen in the numerical results for large m^. 

The data for the critical point coupling, from our different methods of estimating it, separate 
in the parameter range where the latent heat disappears indicating a possible multi-critical 
point or richer structure. For larger values of m? = p we have not been able to resolve the 
nature of the transition. The numerical evidence shows the structure of a 3rd order transition, 
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a behaviour typical of many matrix models, however, the fact that a^; the crossing point of 
the average action curves for different A^, still reliably predicts the transition line, suggests the 
transition is continuous with asymmetric critical fluctuations, consistent with the theoretical 
analysis. 

Our conclusion is that the full model has the qualitative features of the model with m = 
II = (model 5*0 of eq. (13.31) ) and that the effect of the potential is to shift the transition 
temperature and provide a non-trivial background specific heat. The divergence of the specific 
heat arises from the interplay of the Chern-Simons and Yang-Mills terms. Given that this 
competition leads to a divergent specific heat, we expect, sufficiently close to the transition, 
to see the effect of this competition emerge and the specific heat to eventually rise above the 
background provided by the potential and diverge as the critical point. 

The model of emergent geometry described here, though reminiscent of the random matrix 
approach to two dimensional gravity [4] is in fact very different. The manner in which spacetime 
emerges is also different from that envisaged in string pictures where continuous eigenvalue 
distributions [3] or a Liouville mode [23] give rise to extra dimensions. It is closely connected 
to the DO brane scenario described in [12] and the m = version is a dimensionally reduced 
version of a boundary WZNW models in the large k limit [13]. A two matrix model where 
the large limit describes a hemispherical geometry was studied in [32]. It is not difficult 
to invent higher dimensional models with essentially similar phenomenology to that presented 
here (see [19], [21] and [20]). For example any complex projective space CP^ can emerge from 
pure matrix dynamics by choosing similar matrix models with appropriate potentials [19]. 

In summary, we have found an exotic transition in a simple three matrix model. The 
nature of the transition is very different if approached from high or low temperatures. The high 
temperature phase is described by three decoupled random matrices with self interaction so 
there is no background spacetime geometry. As the system cools a geometrical phase condenses 
and at sufficiently low temperatures the system is described by small fluctuations of a U{1) 
gauge field coupled to a massive scalar field. The critical temperature is pushed upwards as the 
scalar field mass is increased. Once the geometrical phase is well established the specific heat 
takes the value 1 with the gauge and scalar fields each contributing 1/2. 

We believe that this scenario gives an appealing picture of how a geometrical phase might 
emerge as the system cools and suggests a very novel scenario for the emergence of geometry in 
the early universe. In such a scenario the temperature can be viewed as an effect of other degrees 
of freedom present in a more realistic model but not directly participating in the transition we 
describe. In the model described in detail above, both the geometry and the fields are emergent 
dynamical concepts as the system cools. Once the geometry is well established the background 
scalar decouples from the rest of the physics and is always massive. If a realistic cosmological 
model the can be found, such decoupled matter should provide a natural candidate for dark 
matter. 
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